Coupling-geometry-induced temperature scales in 
the conductance of Luttinger liquid wires 



p. Wachter^, V. Meden^, K. Schonhammer^ 

^ Institut fiir Theoretische Physik, Universitat Gottingen, D-37077 Gottingen, 
Germany 

^ Institut fiir Theoretische Physik A, RWTH Aachen University and JARA - 
Fundamentals of Future Information Technology, D-52056 Aachen, Germany 

E-mail: waechterOtheorie . physik . uni-goettingen . de 

Abstract. We study electronic transport through a one-dimensional, finite-length 
quantum wire of correlated electrons (Luttinger liquid) coupled at arbitrary position 
via tunnel barriers to two semi-infinite, one-dimensional as well as stripe-like (two- 
dimensional) leads, thereby bringing theory closer towards systems resembling setups 
realized in experiments. In particular, we compute the temperature dependence of the 
linear conductance G of a system without bulk impurities. The appearance of new 
temperature scales introduced by the lengths of overhanging parts of the leads and 
the wire implies a G(T) which is much more complex than the power-law behavior 
described so far for end-coupled wires. Depending on the precise setup the wide 
temperature regime of power-law scaling found in the end-coupled case is broken up 
in up to five fairly narrow regimes interrupted by extended crossover regions. Our 
results can be used to optimize the experimental setups designed for a verification of 
Luttinger liquid power-law scaling. 
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1. Introduction 

Theoretically, it is well established that the low-energy physics of a wide class of one- 
dimensional (Id) metallic electron systems (quantum wires) with two-particle interaction 
is described by the Luttinger Liquid (LL) phenomenology [Ij. One of the characterizing 
properties of LL physics is the power-law scaling of a variety of physical observables 
as function of external parameters. For spin-rot at ionally invariant and spinless systems 
the corresponding exponents can all be expressed in terms of a single parameter, the 
LL-parameter K. On the experimental side the situation is less clear. Much effort has 
been put into measurements aiming at a verification of LL behavior by observing the 
predicted power-law scaling. However, in many of those experiments sources for the 
observed behavior other than LL physics cannot be ruled out in a completely satisfying 
way p]. 

The gap between the status of theory and experiment is partly related to an 
insufficient theoretical modeling of the experimentally available setups. In many 
experiments on quantum wires realized e.g. by single wall carbon nanotubes or 
semiconductor heterostructures, the temperature dependence of the linear conductance 
G was measured [3]. A typical experimental system in which the finite-length LL wire 
is coupled via tunnel barriers to two quasi two-dimensional, stripe-like (Fermi liquid) 
leads is sketched in figure [H In the preparation the precise position of the contact 
regions and its width, as well as the quality of the contacts is still difficult to control. 
In contrast to the experimental geometry the theoretical description is mostly done by 
considering end-contacted wires, in which the lead electrons tunnel into the end of the 
wire and the leads terminate at the contacts. Here we partly bridge the gap between 
the simplicity of the theoretical modeling and the complexity of the experimental setup 
by considering a microscopic model of a quantum wire with Id as well as stripe-like 
(two-dimensional) leads, which arbitrarily couple to the wire. Undoubtedly, in order 
to achieve direct comparison to experimental data true ab-initio simulations of the 
specific experimental setup would be desirable. Yet, the available ab-initio methods 
do not capture LL physics and thus cannot be used in the present context. Therefore, 
microscopic modeling is the method of choice to extend the LL physics originally derived 
within effective field-theories. We show that the length scales set by the overhanging 
parts of the leads and the wire introduce new energy scales, which strongly affect the 
temperature dependence of the conductance. The G{T) curves become much more 
complex than the ones obtained from the simplified modeling. 

Throughout this paper we consider the case of spinless fermions. We expect that 
our results — with the appropriate changes [2] in the equations relating scaling exponents 
to the LL parameter K — directly carry over to the case of spin-1/2 electrons if the 
two-particle interaction is long-ranged (unscreened) in real space (that is, if the gi^± 
term of the g-ology classification [1] can be neglected) and spin-rot at ionally invariant. 
This is generally assumed to be the case in single wall carbon nanotubes [5J. For 
semiconductor heterostructures the range of the interaction can strongly depend on the 
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Figure 1. Sketch of a typical experimental setup. The quantum wire is attached to 
higher dimensional leads, which further extend as indicated by the arrows; note the 
overhanging parts of the leads and the wire. 

screening properties of the surrounding environment (free carriers). In the presence of a 
sizable backscattering component {gi^± > 0) the asymptotic low-energy regime is usually 
reached via a complex crossover behavior [SI [71 [S] . We expect the same to hold for the 
more realistic setups studied here, which would render the behavior of G{T) even more 
complex (and less "universal" ) than described in the present work. 

The most elementary understanding of the temperature dependence of the linear 
conductance of a LL wire can be gained using Fermi's Golden Rule like arguments. 
The local single-particle spectral function p of a translationally invariant, infinite LL 
scales as p ~ (maxjcj, T})"''""^, with u measured relative to the chemical potential and 
«buik = {K + — 2)/2 [2]. For a semi-infinite LL and uj,T vp/x, with x being the 
distance from the boundary and vp denoting the Fermi velocity, the power-law scaling 
of the now x-dependent spectral weight is instead given by the boundary exponent 
ttend = — 1 P, [in]- Beyond the energy scale vp/x the bulk exponent is recovered. 
Fermis Golden Rule then predicts, G{T) ~ T"''""' for tunneling into the bulk of a LL 
while tunneling into the end leads to G{T) ~ T""'^'^. The LL parameter < < 1 (for 
repulsive interactions; if = 1 in the noninteracting case) is model dependent [HE]- Only 
for a few cases (for an example see below) its exact dependence on the model parameters 
such as the strength of the two-particle interaction and the filling of the band is known 
analytically. For small interactions parameterized by the amplitude U, K generically 
goes as K = 1 — U/Uc + 0([t//[/c]^) with the model and parameter dependent scale 
Uc- From this Taylor expansion one can infer that a^nd is linear in U/Uc and for small 
interactions dominates over Ofbuik which is quadratic. 

In experimental setups the particles leave the LL wire at a second contact. To 
obtain the conductance of this type of system one is tempted to simply take the inverse 
of the sum of the two resistances resulting from the tunnel barriers. In the presence of 
inelastic scattering processes this is certainly correct. In the absence of such processes, 
an approximation used in the present paper, it is however less clear if this leads to 
the correct result for the conductance. In [11] and [12] it was shown using scattering 
theory that for low-transmittance contacts and at temperatures sufficiently larger than 
the scale vp/L, with the Fermi velocity vp and the length L of the interacting wire, 
adding the two resistances and taking the invers gives the correct result. We note in 
passing that this can not be extended to the case of more than two impurities [12] . 

The Fermi Golden Rule like considerations on the conductance of a LL wire usually 
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do not take the geometry of the reservoirs into account exphcitely but only use their 
Fermi hquid properties. For the case of end-contacted wires the above results were 
confirmed within a microscopic modeling of the semi- infinite leads (Id noninteracting 
tight-binding chain) and the wire (Id tight-binding chain of spinless fermions with 
nearest-neighbor interaction) [II1[I2]- The linear conductance shows power-law scaling 
j^Qond ^Qj- temperatures Vp/N <^ T <^ B, with the number of lattice sites in the 
interacting wire and the B bandwidth. Other approaches to include the end-contacted 
leads into the model are based on the so-called local Luttinger liquid picture [I3] and 
on the concept of radiative boundary conditions [H], where the wire is modeled by an 
effective field-theory (bosonization) [2]. 

We here extend the functional renormahzation group (fRG) approach used to study 
the linear response transport [H] as well as the finite bias steady-state nonequilibrium 
transport [15j through an end-contacted wire (described by a microscopic model) to 
investigate the role of the contacts and the leads more thoroughly. The reservoirs 
are modeled as two-dimensional tight-binding stripes of variable width (including the 
case of Id leads) and contacted to the Id interacting wire at arbitrary positions. This 
generically leads to overhanging parts of the wire as well as of the leads (see figure 
[1]). We show that they set new energy scales and the temperature range of power-law 
scaling of G{T), vp/N <^T <^ B^ breaks up into up to five different regimes of variable 
size. This leads to a temperature dependence of G{T) which is much richer than the 
simple scaling T"°"<i, which for sufficiently long wires, holds over several decades in T 
within the simplified modehng described above. Our results show that it is difficult to 
find clear evidence for LL behavior of G{T) in the generically used setups. 

To fully cover all relevant temperature regimes and treat systems of experimental 
length, i. e. in the micrometer range, we have developed an algorithm which allows us 
to compute the conductance for systems of up to = 10^ lattice sites. Including the 
stripe-like leads requires a substantial extension of the order algorithm developed for 
Id leads [I6]. 

The rest of the paper is organized as follows. In section [2] the microscopic model 
used to describe the quantum wire and the leads is introduced. In section [3] and 
the appendix it is shown in detail, how to calculate the conductance through the 
system using scattering theory. Section H] is a short survey of the fRG method and 
the generalizations necessary to include stripe-hke leads are described. In section [5] our 
results for the temperature dependence of the conductance are presented and discussed. 
We conclude with a summary in section [61 

2. The Microscopic Model 

We consider spinless fermions on a finite Id lattice with nearest-neighbor hopping t and 
nearest-neighbor two-body interaction f/, coupled to two semi-infinite, stripe-like, and 
noninteracting tight-binding leads. For simplicity all lattice constants are assumed to 
be equal and chosen to be unity. The system is sketched in figure O The Hamiltonian 
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H^irc of the site quantum wire reads (t > 0) 

Hmre -^^kin ~l~ -^int 



with 



N-1 



N-1 



.(t) 



^int = uJ2( - ^) (sVi^.+i - ^ ) , (1) 



where denotes the fermionic annihilation (creation) operator in Wannier states at 
site j G {1, A^}. We here consider the case of half-fiUing of the leads as well as of the 
wire. To assure the latter the density operator Uj = c^jCj in the last line of ([1]) is shifted 
by —1/2. We expect similar results to hold away from half- filling. 
The leads are assumed to be noninteracting 

n,m£A 

where a stands for L (left) and R (right) and the symbol A = {{n^, ny)\nx G {1, Na^^} 
Auy G {1, oo}} denotes the set of lattice points of lead a. The extensions of the leads 
in the x-direction given by iV^ are allowed to differ. The hopping matrix elements ^ 
are assumed to connect nearest neighbors only. 

The Id wire and the stripe-like leads are coupled via hopping terms, described by 

N 

^c"oupl = E E {^f,nC]dH,a + H.C.) . (2) 
j=l neA 

The matrix elements V^"^ will be specified later. 
The full Hamiltonian H of the system reads 

The model for an infinite isolated wire -ffwire with interaction U between all sites can be 
solved exactly using the Bethe ansatz [T7j. The system is a LL for all fillings and 
interactions except for half-filling at \U\ > 2t, where a charge density wave forms 
{U > 2t) or phase separation occurs {U < —2t). From the Bethe ansatz solution, 
the LL parameter K can be computed and reads for half filling and \U\ <2t [T8] 



3. Single particle scattering 

3.1. General relations 

The fRG procedure in the approximation scheme described in section Hlyields a frequency 
independent self-energy of the interacting wire at the end of the fRG flow. This 
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Figure 2. (Color online) Sketch of the Id wire connected to stripe-like leads; the 
crosses depict the sites of the leads, the bars those of the wire; note that the couphng 
between the wire and the leads is located at arbitrary positions leading to overhanging 
parts. 



can be interpreted as a single-particle (scattering) potential and incorporated into a 
renormalized single-particle Hamiltonian. As described in detail in [11] for end-coupled 
wires with Id leads one can compute the conductance from single-particle scattering 
theory [19j using a generalized Landauer-Biittiker approach |20] . 
We therefore consider the single-particle Hamiltonian 

h = h„ire + hi^nd + ^fead + ^coupl + ^coupl • (5) 

where /^wire is the effective single-particle Hamiltonian for the wire at the end of the fRG 
flow (which includes the generated scattering potential) and the /I's are single-particle 
versions of the terms introduced in the last section. For the states in the single-particle 
Hilbert space we use standard Dirac notation. In the scattering problem to be solved, the 
first three terms of h are considered to be the unperturbed part ho, while the couplings 
to the leads /icoupi + ^mupi considered as the perturbation hi. We closely follow the 
steps for the corresponding scattering problem with end-contacted. Id leads [TT] . 

The main difference of working with higher dimensional leads is the occurrence of 
transverse scattering channels. The orthonormal single-particle eigenstates \k,l,a) of 
the isolated semi-infinite leads are standing waves with the wavenumber k G [0, vr] in 
the y-direction and / = 1, N^^a labeling the transverse modes 

{n,a\k,l,a) = \l — sin{kny) 




In 

■ sm I — 



= {ny\k){n^\ir = ''{n\kr. 
The corresponding eigenenergies are 

e,(fc, /) = - cos(/c) - 2tl cos (j;^^^^^ 

= ea{k) + ea^i = ea{k) , (6) 
with the nearest-neighbor hopping matrix elements t% and 
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We now consider the scattering of a particle from the left to the right lead (or vice 
versa) coupled via hi to the wire. The outgoing scattering states \k,l,a+) = \k+)°- 
follow from the Lippmann-Schwinger equation [21] 

\k+)-=\k)- + g^ea{k)+iO)hi\kr , (7) 



where Q{z) = {z — h) ^ is the resolvent of the full single-particle Hamiltonian. 
If we define a as the "complement" of a this yields 

-^{m\k,+r = J2"{m\g{ea{k)+tO)\j)V;/{n\kr (8) 

Using g = Go + GohiQ with Qo{z) = {z — ho)^^ we obtain 

-^{m\k,+y = 5^^(m|6;o|mrv;i',,'(/|6^|j)^^"(r?|fc)^ (9) 



3d' 

rh' ,ri 



Only elastic scattering k,l,a ^ k', I', a occurs with ea{k, I) = ea{k', I'). This is enforced 
by the free resolvent matrix elements "'{m'\Qo{ea{k, I) +iO)\rn)°-. It is easier to see this in 
a mixed representation, local in the ^/-direction and using the transverse channel number 
/' as the additional quantum number 

'^{my,l'\k,+y = 

3,3' 



TYh'y 



x^{n\kr)^=. (10) 



The second equality holds for niy larger than the involved in the coupling of the lead 
and the quantum wire. To derive this relation we used the result for the resolvent of a 
semi-infinite Id chain [n = 1, 2, 00) with nearest-neighbor hopping t 

_ 2sin[A;(e)n] .^^^^^ 

v[k{e)] ' ^''^ 

where the second line holds for m > n. For |e| < 2t the wavenumber k{e) G [0, vr] is 
the solution of — 2t cos[A;(e)] = e and the velocity v{k) is given by 2tsin(/c). In |Tl] only 
the result for n = 1 was presented. For |e| > 2t the resolvent matrix element decays 
exponentially with the distance \m — n\ and are thus "closed" channels in ( |TOl) . This 
happens if \ea{k,l) — ea,i'\ < 2|t^|, and k' G (0, vr) is determined by energy conservation 
eaik') = ta{k) + ea,i - ea,i'. The velocity in ([10]) is given by v^{k') = 2t^sm{k'). The 
prefactor of the outgoing plane wave in (JTOl) gives the transmission probability for the 
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scattering k,l,a ^ k', I', a 



rh.fi 



(12) 



Taking the sum over / compatible with ea{k) = e and over the /' values of the 
corresponding "open" channels we define pO] 



Ue) = J2^^:^\t{eJ,l')\\ (13) 



{1,1'} 

The unitarity of the S'-matrix guarantees Ta{e) = Ta{e) = T{e). This quantity 
determines the conductance G{T) in the generalized Landauer-Biittiker formula [20] 



G{T) = ^ /de (-^ ] T{e) , (14) 



h J \ de, 

where / denotes the Fermi function and /h is the unit of the quantized conductance. 
Note that for our interacting problem the fiowing self-energy and thus /i„ire as well as 
the effective transmission probability T{e) become T and N dependent (for details, see 
below) . 

The Landauer-Biittiker formula f|T^ holds in the abscence of inelastic processes 
due to the two-particle interaction (that is for vanishing imaginary part of the self- 
energy) [19]. Such processes do not appear in our approximate treatment of the 
interaction and using f|T^ does thus not present an additional approximation. 

For the calculation of the transmission probability one has to compute the matrix 
elements {j'\Q\j) of the full resolvent. This can be reduced to the inversion of a x A^- 
matrix, using Feshbach projection [TT] 

pg{z)p = 

[zP - PhP - PhQ{zQ - QhQ)-^QhP] , (15) 

where P and Q = 1 — P are projection operators. Here we use as P the projection 
operator Pwire = 12f=i onto the wire and Q the one on the leads. 

The special geometry (mimicking a generic experimental setup) used in section 
[5] is shown in figure [3l Each site of the wire lying on top of a site of the lead (as 
shown in the figure) is coupled to this site via a hopping term fr^. The length scales 
appearing in the system are also defined in the figure. Note that the overall length 

The projected inverse Green function takes the form shown in figure lAll of the 
appendix. We need the tridiagonal elements of this Green function to calculate the 
right-hand side of the fiow equations (ITSl) and (fT9l) . see section HI An algorithm to 
compute these elements is presented in the appendix. Furthermore, to calculate the 
transmission amplitude ( |T2i) we need the j, j'-elements of the Green function with 
j' G {Ni + 1,..,A'2} and j G {N^, + 1,..,A^4}, i. e. the elements connecting the left 



Coupling-geometry-induced temperature scales in the conductance of LL wires 



9 



X X X X 1 „ 

XXXIfvrL XXX X Nl 

X X X J lead X X X X J "^'"^ 

■ I I I I 



No^g, XXX N^y X X X X N„„ 

XXX X X X X 

XXX X X X X 



Figure 3. (Color online) Sketch of the system considered. The overlapping sites of 
wire and leads are individually coupled via a hopping term. The length scales are 
indicated. 



and right contact region. An algorithm to efficiently compute these elements is also 
briefly sketched in the appendix [22]. 

For simplicity, we assume the hoppings in the leads and in the wire to be equal and 
define this as our unit of energy, i. e. t'^ = ty = t = 1. 

3.2. One- dimensional leads 

The appearance of new temperature scales is most easily seen in the case of one- 
dimensional leads, i.e. A^^^^^ = = 1. For ^ ^lead ^^^1 leads weakly coupled 
to the wire, i.e. [t^'^)^ <^ 1 the transmission probability T(e) consists of A^ narrow 
resonances near the eigenvalues ea of the wire without contacts. The integrated weight 
Wa of the Lorentzian resonance near is [TI] 

-1 



1 1 



(16) 



Here the partial widths are given by 

2 



= {t'^y 



{N + l)sm[k{e^)] 

with A^" the site of the wire to which lead a couples. In the temperature regime 
Vp/N <^T <^ B the Landauer Biittiker formula simplifies to 



( df 



G(r)«^»„(-^j (17) 

We return to this expression when discussing our numerical results. 



4. The Functional Renormalization Group 

In this section we present a short outline of the fRG at finite temperature T, which we 
use to treat the interacting system. General aspects of the fRG can be found in [23], [21], 
and [25j, while its application to inhomogeneous LLs is described in [8], [11] and |16j . 
By comparison to results for small systems obtained by density-matrix renormalization 
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group and to exact results from bosonization and Bethe ansatz for the asymptotic 
behavior, it was shown that the fRG in the truncation scheme used here captures the 
relevant physics not only in the asymptotic low-energy regime but also on finite energy 
scales [TT1[T6]. Towards the end of this section, we describe the extensions of the method 
necessary when dealing with stripe-like leads. 

Starting point of our implementation of the fRG is the introduction of an infrared 
energy cutoff A in the Matsubara frequencies of the noninteracting propagator. This 
leads to a A-dependence of the generating functional of the one-particle irreducible Green 
functions (vertex functions). Taking the derivative of this functional with respect to A 
leads to an exact, but infinite, hierarchy of coupled first order differential equations for 
the vertex functions. We truncate this hierarchy by neglecting the m-particle vertices 
with m > 3. In order to reduce the numerical effort, we project the 2-particle vertex 
function onto the Fermi points and parameterize it by a nearest-neighbor interaction 
of amplitude t/^ [HlIIS]. Finally, we neglect the feedback of the 1-particle vertex (the 
self-energy) onto the flow of the 2-particle vertex function, that is the flowing effective 
interaction [16j. The resulting set of equations is integrated from the initial cutoff 
A = oo down to A = 0, at which the original, cutoff-free problem is recovered. For end- 
contacted, interacting wires and weak to intermediate interactions this approximate 
approach was shown to provide reliable results by comparison with exact analytical and 
numerical results [HI [HI [IE] . 

Within our approximation the self-energy S becomes real, frequency-independent 
and tridiagonal in position space. It acquires a dependence on the temperature and size 
of the interacting wire and can thus be interpreted as a (T- and iV-dependent) single- 
particle scattering potential. This leads to an effective single-particle Hamiltonian with 
renormalized onsite energies and hopping matrix elements. 

As mentioned in the introduction to leading order the scaling exponent of the 
tunneling conductance into an end-coupled LL is linear in U. In contrast, tunneling 
into a translationally invariant (bulk) LL is characterized by an exponent which for 
small U is quadratic p]. First order scaling exponents were also found in other setups 
of inhomogeneous LLs (e.g. LLs with a single impurity [9l jTTl |T6] and Y-junctions of 
LLs; see e.g. Ref. [26] ). In all cases in which predictions for the scaling exponents of 
inhomogeneous LLs from alternative approaches (e.g. effective field-theories) exist those 
were reproduced to linear order using our truncated fRG method [HI [161 [26]. However, 
our truncation procedure does not capture bulk LL physics, as terms of order f/^ are 
only partly kept in the flow equation for the self-energy. In particular, the terms of order 
U"^ In (uj) obtained in a perturbative calculation of the self-energy of a translationally 
invariant LL (evaluated at the Fermi momentum kp), leading to the bulk LL power 
laws when resummed by an RG method, are neglected in our approach. We thus expect 
that, in situations which are dominated by bulk LL physics (for examples see below) 
the scaling exponent vanishes within our approximation. 
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The set of flow equations to be solved reads 

— S^. = - - 
dA ^'^ 71 

AyA = i 

dA TT 



^S^. = -^Ef^i^>^f^^[^k.+.(^A)] , (18) 
±1 

4K^tm = zU'.^i^-Kmm^ (19) 



f/^ = 7 ^ , (20) 



v^i+A^y 27r 

where we specialized the flow equation for the effective interaction to the case of half- 
filling. The implementation of finite temperatures requires at each step of the flow the 
replacement of the flow parameter A by the nearest Matsubara frequency [8]. 

In order to compute the right-hand side of the flow equations for one has to 
compute the tridiagonal elements of the matrix Q^{iA) inverting the given N x N- 
matrix Q^{iA)^^. To cover all relevant energy scales and to treat systems of lengths in 
the micrometer range we need an algorithm, which allows us to calculate these matrix 
elements more efficiently than the standard order A^^ methods [27]. In the case of Id 
leads, Q^{iA)~^ itself is tridiagonal and the tridiagonal elements can be computed in 
order N time [161 l2H]- For stripe-like leads Q^{iA)~^ consists of tridiagonal blocks, 
connected on the sub- and superdiagonal to full blocks, one for each contact region. 
The explicit structure is depicted in figure lAll of the appendix. There we present an 
algorithm, in which the problem of inverting the full A^ x A^-matrix is split up into the 
multiple inversion of each isolated block. Thus, we can readily use the order A^ algorithm 
for the tridiagonal parts and some standard algorithm [27] for the full matrices. The 
bottleneck of the method is therefore the width of the stripe-like leads, which determines 
the size of the full matrices. 



5. Results 



5.1. One- dimensional leads arbitrarily coupled to a noninteracting wire 

We first consider the conductance G{T) through a clean quantum wire with two semi- 
infinite Id leads (stripe-like lead with a width of one lattice site) coupled at arbitrary 
position of the leads and the wire. This introduces four additional length scales, namely 
^leaif ^^"^ NoJef, apart from the overall length A^ of the wire (see figure [3]). They strongly 
effect the temperature dependence of the conductance. 

First we shortly comment on even-odd-effects arising from the fact that the coupling 
sites in the wire as well as those in the leads can each be even or odd and so can be the 
overall length of the chain. Within the notation introduced in section 13.11 the coupling 
sites in the wire are A'^gj. + 1 and A't^er + those in the leads A'lead + 1 N^^^ + 1. 
For T — 0, the derivative of the Fermi function in f|T^ becomes a 5-function and the 
integral is given by the effective transmission evaluated at the Fermi energy. This shows 
even-odd-effects. For T — the conductance reaches a nonvanishing value only if the 
coupling sites are all odd. For even sites the local spectral density at the Fermi energy 
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Figure 4. (Color online) Conductance G as function of temperature T of a 
noninteracting wire of length N — Iff* + 1 and coupling t"^^ = t'^^ — 0.25. Solid 
line: end-contacted wire, i.e. NoJJx = (setup shown as the right inset) obtained from 
([Ti)) . Circles: TV^^j, = 20, TV^^j, = 286 (setup shown as the left inset) obtained from 
([T^ . Dashed line: the same parameters as for the circles but using the approximate 
analytical result ([TT]). The vertical lines indicate the crossover scales resulting from 
the overhanging parts. 

vanishes, leading to G{T = 0) = 0. If, for odd coupling sites and t^^ = t'^^, the overall 
length is taken to be odd, G(T — 0) — e'^/h holds. For temperatures of the order of 
vf/N and higher these even-odd effects vanish. We mainly consider systems with odd 
and coupling at odd sites. This implies G{T = 0) = e'^/h for f^^ = f^^, which also 
holds for the interacting system. 

In figure HJ the conductance G(T) of a noninteracting wire coupled at the ends is 
compared to the conductance of a system coupled in the bulk. The leads are weakly 
connected with t^'^ = f^^ = 0.25 and terminate at the contacts (no overhanging parts, 
N^L^^ = iVif^^ = 0). For the bulk-coupled case we chose N^^^^ = 20, N^^^. = 286 (circles). 
The conductance of the end-coupled wire (solid line) starts at e'^/h for T = 0, crosses 
over into a plateau at a scale T ^ vp/N (left vertical dashed line) and finally falls off 
as T"^ for T ^ i? (band effect). For the bulk-coupled case the plateau region breaks 
up into three clearly distinguishable regimes with the two additional crossover scales 
vp/N^^^ and Vp/N^^^^ indicated by the other two vertical dashed lines. In the first 
region with vp/N <^ T <^ vp/N^^^ a plateau with the same conductance as for the 
end-coupled chain appears. At the two additional scales the conductance crosses over to 
new plateaus with lower conductance values and finally falls off for T 3> -B. The dashed 
line shows the conductance of the bulk-coupled system using the approximate analytical 
expression (|T7|) for comparison. For temperatures in the range vp/N<^T<^Bwe find 
excellent agreement with the exact result obtained from ( |T4l) . 

Note that the overhanging parts of the wire and of the leads enter symmetrically in 
ffTTll . Thus, choosing A^i^^^ = 20, N^^^^ = 286 and A^^^^, = A^^^^, = 0, gives an identical 
conductance as that shown by the circles and dashed line in figure |H Obviously this 
will no longer hold in the presence of interactions in the wire. 
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Figure 5. (Color online) Conductance G of a noninteracting wire of length A'' = 10^+1 
and coupling t'^'^ = = 0.25. Solid line: end-coupled wire, i.e. N^J^^ = NoJcf ~ 
obtained from the analytical result (fT7|) . Dashed line: N^^^^ ~ 10^,iV^gj. ~ 2 • 
W'.N^^^ = 2086 and N^^^ = 86 (setup shown as the inset) obtained fr om the 
analytical result (fT7|) . The vertical lines indicate the crossover scales. 

All four new temperature scales, i.e. vpjNyJJ^ and vp/NoJ^ can be resolved using 
the analytical expression ( fTTI) for a very large system (that is for a very small lower 
bound vp/N). This is shown as the dashed line in figure [51 The solid line shows the 
conductance for a system coupled at the ends for comparison. 

After identifying the relevance of the additional energy scales due to the overhanging 
parts of the wire and the leads for noninteracting wires we turn to the interacting case. 

5.2. One-dimensional leads arbitrarily coupled to an interacting wire 

As described in section H] the approximate fRG flow leads to a nontrivial static self- 
energy for the interacting wire which is tridiagonal in the site indices. At half-filling 
and for tunnel barriers as used here the diagonal part vanishes and the effective potential 
is given by a nonlocal term, the modulated hopping j+i. If the leads are coupled to 
the ends of the wire (with and without overhanging parts of the leads) , the two contact 
barriers are the source of (nonlocal) potentials which oscillate as (— 1)-' (for half-filling) 
and decay as where j denotes the distance from the ends (at j = I and j = N, 
respectively) . The prefactor of the decay is universal in the sense that it does not depend 
on the strength of the tunnel barriers t'^'^ and f^^ and is the same even for t'^'^ = f^^ = 0, 
that is a decoupled wire with open boundaries. The power-law decay is cut off at a scale 
jr oc 1/T set by the temperature beyond which S^j+i decays exponentially. Scattering 
off this effective potential leads to the power-law scaling of the local spectral weight and 
the conductance discussed in the introduction with the exponents determined by the 
universal prefactor jJl HI] . 

A typical off-diagonal component of the self-energy Sjj+i of a bulk-coupled 
interacting wire is shown in figure El From the ends of the wire the {—ly/j oscillations 
with universal prefactor emerge and become exponentially damped at a distance 
jr oc 1/T (indicated by the vertical dashed line). For clarity only one end is shown in 
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Figure 6. Off-diagonal self-energy at the end of the fRG flow for an interacting 
wire with U = 0.5 and couplings t^^ = = 0.25 at T = 4 • 10"^. The wire of 
length TV = 10"* -I- 1 is coupled to the end of the leads (i.e. N^^^^ = N^^^^ = 0) with 
^over = 1086, N^^j. — 286. The part shown in the main plot is indicated by the dashed 
box in the inset. Long-ranged oscillations emerge from the boundaries and from the 
coupling site. The scale Jt beyond which the oscillations are damped exponentially is 
indicated by the vertical dashed lines. 



the main plot. The couphng to the lead (at site 1087) introduces similar oscillations but 
with a nonuniversal prefactor. In contrast to the situation found close to an end-contact 
these oscillations do not lead to a power-law suppression of the spectral weight close to 
the coupling site resulting from the local inhomogeneity. Furthermore, as discussed 
in section H] our truncated fRG procedure does not capture bulk LL power laws with 
exponents of order [/^. Thus, within our approximation the spectral weight close to the 
coupling site in the bulk shows no power-law behavior at all. 

Figures [7] to [9] show the conductance for several parameter sets and temperatures 
larger than vp/N, i.e. temperatures in the regime in which in the noninteracting case the 
plateaus emerge. The length of the wire is = 10^+1 sites which corresponds to roughly 
a micrometer taking typical lattice constants. The vertical dashed lines terminating at 
the different curves indicate the new energy scales for the corresponding parameter set. 
The curves for the corresponding end-coupled systems without any overhanging parts 
are included in the figures for comparison (dashed lines). The effective exponent of 
GiT), computed as the logarithmic derivative of the conductance, is shown as the right 
inset in each figure. The left inset shows the setups of interest. 

In figure [7] we focus on wires with N^^^j. = N^^^ = 0. For all three temperature 
regimes, 

vp/mi^^iNl^eid) < ^ < VF/mm{Nl^JJ) and 
VF/mm{Ni/J) « T « B, 

the usual power law G{T) ~ T"™'' is found. Note that for the system with N^^^ = 20, 
-^lead — 1086 (dotted line) the first temperature regime vf/N -C T <^ "^FZ-^ifad too 
small for a plateau in the effective exponent to develop. However, by reducing N^^^ one 
can extend this regime, such that the effective exponent tends towards Oend- 

Already this simple case with no overhanging parts of the wire (that is the particles 
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Figure 7. (Color online) Main plot: Conductance G of an interacting wire with 
\J = 0.5 of length iV — 10'' + 1 as function of the temperature T. The coupling to the 
leads is located at the end of the wire {N^^^j. = N^^^. = 0) and has a small amplitude 
^CJ, ^ ^cz, ^ 0.25. Dashed line: TVi^/j^ = 0. Dashed-dotted line: N^^^^ = 0, 7V,^^d = HO- 
Dotted Hue: N^^^ ~ 20, A^i^^d ~ 1086. The vertical dashed lines terminating at 
the different curves indicate the corresponding crossover scales. Left inset: Setups 
studied. Right inset: Effective exponents. The solid horizontal line indicates the fRG 
approximation for aond as obtained in [llj . 



tunnel into the end of the interacting wire) demonstrates our central result. The power 
law, which for end-contacted wires of experimental length without any overhanging parts 
of the leads was found to hold over several decades, is split up into significantly smaller 
temperature regimes (of roughly one decade) with power-law scaling interrupted by 
extended crossover regimes. Thus the three plateaus of different conductance obtained in 
the noninteracting case develop into power laws with the same exponent in the presence 
of two-particle interactions. For each temperature regime the appearance of the latter 
can be traced back to the spatial dependence of the self-energy Sjj+i. Over the whole 
temperature range Vp/N <^ T <^ B, when moving from the left to the right contact the 
fermions have to pass the long ranged oscillations with universal amplitude originating 
from the ends of the wires. For vp/N ^ T these are well separated, such that the 
simple arguments presented previously hold, implying that G{T) ~ 

The role of the length scales resulting from overhanging parts is further exemplified 
in figure [8] for systems with N^^^ = N^^^ = and N^^^ = but varying N^^^^. Beginning 
with a purely end-coupled system (iVj^gj, = 0; dashed line) the left coupling site is moved 
into the bulk of the wire. For temperatures T -C Vp/N^ver particles pass both 
oscillatory potentials originating from the ends of the wire and the situation of fermions 
tunneling in and out of the two ends is effectively recovered. This leads to power-law 
scaling of the conductance with exponent aend- With increasing N^^^^. the crossover scale 
'Vp/N^ver decreases and the effective exponent reaches aend inside a smaller and smaller 
temperature regime (see the inset). For the system with the largest iVj^g^. (iVj^gj. = 1111; 
dashed-dotted line in figure [HI), the interval [vf/N,vf/N^^^^] becomes too small for the 
logarithmic derivative to reach its asymptotic value (see the inset). For temperatures 
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T > vf/N^^^^ the oscillat ions originating at the left end of the wire are exponentially 
damped at the position of the left contact and to the right of it. On their direct 
path from the left to the right contact the fermions thus only experience the potential 
originating from the right end of the wire. As discussed in the introduction one such 
potential is sufficient for the spectral weight close to its origin to scale as 0;°™'*. One 
is thus tempted to conclude that G(T) follows a power law with exponent ttend also for 
'^p/N^ver ^ T <^ B. Figure [8] shows that this is not the case. Instead, the effective 
exponent seems to approach an asymptotic value which is significantly smaller than 
«end (see the dashed-dotted line). The naive expectation ignores the complexity of the 
scattering problem. For example, the fermions can also move to the left of the left 
contact before reaching the right lead. Below we further characterize the behavior in 
this regime. 
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Figure 8. (Color online) As in figure [7] but with N^^^^ = TVj^^^ = 0, N^^^ = and 
for different N^^ci- Dashed line: N^^^^. — 0. Dotted line: N^^^^ — 11. Dashed-double- 
dotted hue: N^,,,, = 111. Dashed-dotted Une: N^„,, = 1111. 



As a third example in figure [9] we show results obtained for N^^^ = N^^^ = 
(no overhanging parts of the leads) and different positions of the coupling to the wire 
NoJer'- In the regime Vp/N <^ T <^ vp /'s\ax.{NoJJ^)^ the conductance follows the power 
law G{T) ~ T"™<i since the oscillations of '^jj+i from the two boundaries reach beyond 
the positions of the contacts. For N^^^ = 1086 (dotted line) this regime becomes too 
small for the power law to develop. For vp/max{NoJe?') ^ T ^ VF/viim{NoJer'), the 
oscillations from the boundary, originating at the site of the longer overhanging part, 
are already cut off as the exponential damping sets in at Jt oc 1/T, i.e. Jt <^ max(A''OTef ). 
In this regime the curves seem to show power-law behavior with an exponent smaller 
than Oend similar to the behavior discussed in connection with figure [HI In the third 
regime vp/mm{NoJef) <C T <^ 5, the oscillations from the boundaries do not reach 
the region between the left and right contact (due to the exponential damping) and 
no power-law behavior can be found (see the conductance plateau of the dotted and 
dashed-dotted line in the main part of figure [9]) . This is in accordance with the fact 
that the oscillations of ^jj+i originating from the bulk contacts do not have the same 
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effect on the spectral weight as the ones coming from the boundaries of the wire (no 
power-law suppression). We expect the conductance in this temperature regime to show 
power-law scaling with the bulk LL exponent Obuik ~ U"^ not captured by our truncated 
fRG procedure. 
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Figure 9. (Color online) The same as in figure [7] but with N^^^ — N^'^^ — 
and coupling into the bulk. Dashed line: Novif — 0. Dashed-dotted line: N^^^^ = 
20, N^^^ = 286 Dotted line: A^i;„. = 20, N^^^ = 1086. 



To obtain a better understanding of the apparent second exponent smaller than 
ctend found above, figure [TOl shows the logarithmic derivative of G(T) in the temperature 
regime vp/msixN^^^j. <^ T <C i? for a system with N^^^j. = 1110 and N^^^. = 0. Curves 
for different interaction strengths U are plotted. For U < 0.5 (double-dashed-dotted hne 
and below) one is tempted to conclude that G scales like T"™<^/^. The lower horizontal 
line shows aend/2 for U = 0.5 numerically calculated within the fRG-implementation 
used here [H]. However, for U > 0.5 the plateau of the logarithmic derivative becomes 
strongly tilted, such that no definite statement about any power law can be made 
(compare the solid horizontal line showing aend/2 for U = 1 with the dashed curve). 
Attempts to treat the underlying scattering problem, e.g. with phase- averaging [I2] 
or methods developed for Y-junctions [26] have not yet led to a consistent picture. 
Furthermore, in this temperature regime G{T) might be altered if bulk LL behavior is 
properly introduced. 

The conductance through "mixed" systems, that is those with overhanging parts 
of the quantum wire as well as overhanging parts of the leads, shows an even richer 
temperature dependence. As in the noninteracting case (see figure [5]) up to five different 
temperature regimes and the corresponding crossover regimes might appear in 
[vf/N,B]. The behavior sufficiently away from the characteristic energy scales set by 
the lengths of the overhanging parts can be understood as explained above. 

These considerations show that for a generic experimental system with four 
overhanging parts all having different lengths and not all of them being several orders 
of magnitude smaller than the total length of the wire it will be very difficult (if not 
impossible) to observe clear indications of a single power law in G{T). We next show 
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that considering stripe-like leads, as being realized in the experiments, does not improve 
this situation. 

5.3. Stripe-like leads arbitrarily coupled to a quantum wire 

We here extend our analysis to stripe-like leads, as realized in many experiments. For 
this setup the parameter space grows linearly with increasing the number of coupling 
sites, e.g. for five sites at each contact, we can freely choose ten couplings t'j^^^. To 
prevent a proliferation of parameters we therefore focus on equal couplings, tj^^"^ = t . 

Figure [TT] shows the conductance of a noninteracting quantum wire with A'mid = 
963, N^^^^ = N^^^^ = 18, A^i^^d = N^^^^ = and for varying width of the leads A^'c^o^n^. 
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Figure 10. (Color online) System as in ligure[8]with A^^or = 1110. Effective exponents 
in the regime vp/N^^^^ <C T ^ _B for different interaction strengths U. From bottom to 
top: U — 0.25; 0.5; 0.75; 1.0; 1.5. The horizontal lines show the fRG approximation for 
<^end{U = 0.5) /2 and acnd(C^ — l-0)/2 as obtained in [TT]. The vertical line indicates 
the crossover scale. 
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Figure 11. (Color online) Conductance G as function of temperature T of a 
noninteracting wire with stripe-like leads of different width. All wire-lead couplings 
are set equal: f"^ = f"^ = 0.25. The setup is shown in the inset: TV^^,,. = N^^^ = 18, 
-^load ^ ^load ^ ^ ^"'^'^ N^ud — 963. From bottom to top, the width of the contact 
region is N^^^ = = 1, 3, 5, 7, 9, 11, 15, 21. The vertical lines indicate the crossover 
scales. 
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Figure 12. (Color online) As in figure [TT] but with length A^mid = 10'* + 1 and 
width N^^^ = N^^ = 5 of contact region. Dashed-dotted Una: N^^^ = 0, TV)^^^ = 
and iV^^, = 20, N^^^ = 286 (setup shown in left inset). Dashed hne: N^^^^ = 20, 
-^load — 286 and N^^^j. = N^^^j. — (setup shown in right inset). SoHd line: 
N,^ , = N,^^ , — and = — Dotted line- — — 1 

-^lead = -^lead = ^ and N^^^^ = A^^^^jj. = for comparison. The vertical lines indicate 
the crossover scales. 



The vertical lines indicate the temperature scales vp/Nover and vp/N, which turned out 
to be relevant for Id leads. They seem to be of equal importance for stripe-like leads 
(see below). The overall shape of the curves is unaltered by increasing the width of the 
leads. For large numbers of lead channels G approaches the unitary conductance e'^/h 
of a single channel wire. This trend can be counteracted by reducing the strength of 
the lead-wire coupling. 

In the case of Id leads we identified two pronounced even-odd-effects. The 
dependence on the total length of the wire is also found for stripe-like leads. The 
vanishing local spectral density at even sites of the wire gives rise to the second even- 
odd-effect. For stripe-like leads one coupling site of the wire has to be odd and this 
effect vanishes as the leads become broader. 

Next we further investigate the relevance of the energy scales identified for Id leads 
starting with the noninteracting case. Figure [12] shows the conductance for similar 
configurations as in figure la but with NcJn^ = 5. For comparison, the results for the 
end-coupled system with N,^,^ = N,'^.^^ = 0, A^^^,, = A^^^^, = 0, and Ar^^„, = A^^^^ = 5 are 
included. Overall, the conductance is considerably higher than in figure H] as expected 
from the results shown in figure [TTl At the temperature scales Vp/N, vf/N^^^^.^^^^, and 
'^F I ■^aver/iea.d^ couductauce subscqueutly crosses over to lower plateaus. However, 
the complete equivalence of -/V^gj,^jgj^^ and Nj^^^j^^^^ is lost even for the non-interacting 
case as the geometric equivalence is no longer present. 

Figure [13] shows a typical off-diagonal component of the self-energy for the case 
of stripe-like leads coupled to the bulk of an interacting wire. For comparison to the 
case of Id leads the parameters are chosen the same as in figure [6] The oscillations 
emerging from the boundaries are expectedly unaffected by the width of the leads. The 
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Figure 13. (Color online) Off-diagonal self-energy at the end of the fRG flow for an 
interacting wire with U ~ 0.5 at T = 4 • 10"'^. All wire-lead couplings are set equal: 
^cn ^ ^ 0.25. The setup is shown in the inset: N^^^ = 5, N^^ = 1, iV^id = W^ + l, 
^load = ^icad = and N^^^^ = 1086, iV„^„, = 286. The part shown in the upper plot 
is indicated by the dashed box in the inset. Long-ranged oscillations emerge from the 
boundaries and from all coupling sites. The lower plot is a magnification of the left 
contact region. 

contact region itself (of size A^^^^^; lower panel in the figure) contains a superposition of 
oscillations emerging from all coupling sites. The symmetric decay of these oscillations 
results from the fact that we chose all couplings to be equal. The oscillations originating 
from the left contact to the left and the right contact to the right however still fall off 
as the inverse distance up to a scale Jt oc 1/T beyond which they decay exponentially. 
As for Id leads these oscillations do not induce a power-law suppression of the spectral 
density within our fRG approximation. Because of these similarities the temperature 
dependence of G for an interacting wire with stripe-like leads can be traced back to the 
case of Id leads. 

In figure [14] we focus on end-coupled systems (no overhanging parts) with f^^ = 
= 0.1 (weak coupling). The only energy scales emerging are Vp/N and B. For 
vp/N -C T <^ 5, the conductance scales as ~ T°cnd ^\^q closeness of the results 
for NcJJ^ = 5 and 6 shows the minor importance of the even-odd effect. Figure [T5l 
shows systems with tunneling into the end of the wire but with overhanging leads. For 
comparison the end-coupled case without overhanging leads is shown as the dashed line. 
As in the case of Id leads the universally decaying oscillations from the ends of the 
system induce power laws G{T) ~ T"<="'^ in the three temperature regimes introduced 
by the overhanging parts. For a wire of roughly a micrometer in length it is difficult 
to clearly resolve these power laws due to the extended crossover regions and the upper 
(bandwidth B) and lower (cutoff scale vp/N) bounds of any power-law scaling. We 
thus choose N^^^ = N^^^ in order to reduce the number of regimes to two. Even in 
this case the two power laws cannot be resolved for a single parameter set (compare the 
dashed-dotted and dotted lines). 
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Figure 14. (Color online) Main plot: Conductance G of an interacting wire with 
interaction strength U = 0.5 as function of the temperature T coupled to stripe-like 
leads of different width. The couplings to the leads are set equal: t'^'^ = =0.1. The 
setup is shown in the left inset: N^id = lO'* + 1, A^i^^d = ^icad = 0= ^o^vcr = ^o^r ^ 
and N^^^ = = 1;6;5;11 (from bottom to top as indicated in the plot). The 

vertical dashed lines indicate the crossover scales. Right inset: Effective exponents. 
The solid horizontal line indicates the fRG approximation for aond as obtained in [TT] . 
Line styles as in main part. 
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Figure 15. (Color online) As in figure [T4l but with different iVj^^^ — N^'^^. Dashed 
line: N^^^^ = N^^^^^ = 0. Dotted line: N^^^^ = N^^^^ = 100. Dashed-dotted line: 
^^icad = <ad = 332. 



The case of leads coupled to the bulk of the interacting wire and terminating at 
the contacts is shown in figure [HI As anticipated from the noninteracting case (figure 
[T2l) . the temperature scales introduced by the additional length scales iV^^j. and N^^^ 
are still present. In all three temperature regimes 

vf/N < T < VF/max{N^J^f), 
VF/max{N^J^f) < T < VF/mm{N^JJ^) and 
VF/mm{N^i^) <^T<^B, 

the conductance G{T) is governed by the same behavior as in the case of Id leads. In the 
first regime the oscillations from the boundaries strongly superpose the coupling sites, 
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Figure 16. (Color online) As in figure [14] but with different N^^^^ = N^^^. Dashed 
line: iV^^^ = N^^^ = 0. Dashed-dotted line: N^^^^. = N^^^. = 1086. Dash-dot-dotted 
line: TV.^.^, = 1110, 7V„^,, = 0. 



thus we find G(T) ~ T"^"'^. Note that to better resolve the second and third regime, this 
first regime is chosen such small in the curves shown that the asymptotic value of the 
exponent is not reached. The second regime is exemplified by the dashed-double-dotted 
curve. After crossing over at scale vp/Novcv, the conductance seems to follow a power 
law. However, as explained in detail in the section on Id leads, this can not definitely 
be confirmed and the behavior might be altered in the presence of the bulk LL power 
law. In the third regime, the oscillations in the self-energy from the boundaries are cut 
off by the exponential fall-off at scale oc 1/T before reaching the contact region, thus, no 
power-law behavior can be found within the truncated fRG (figure 16, dashed-dotted 
line). We expect the conductance in this regime to show scaling with the bulk LL 
exponent if computed with an improved approximation which properly describes also 
bulk LL behavior. 

The conductance of "mixed systems" with overhanging parts of the leads and the 
wire again shows an even richer temperature dependence which can be understood by 
combining the cases studied above. 

6. Summary 

In the present study we were aiming at a more realistic modeling of transport through 
interacting Id quantum wires showing LL physics by considering stripe-like leads of 
variable width coupled to the wire at arbitrary positions via tunnel barriers. We showed 
that the overhanging parts of the leads and the wire, which generically appear in 
experimental setups designed for a verification of LL power-law scaling, strongly modify 
the temperature dependence of the linear conductance. For end-contacted wires of 
experimental length with leads terminating at the contacts the range of possible power- 
law scaling (with exponent Oend) bound by Vp/N and B extends over roughly four 
decades. In the presence of up to four additional energy scales induced by the lengths 
of the four overhanging parts and the corresponding extended crossover regimes the 
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regions of possible scaling become so small that the power laws are not visible. This 
shows the difficulties to clearly confirm power-law behavior in transport experiments 
with LL wires. 

One might speculate that the situation becomes even worse if inelastic two-particle 
scattering processes neglected in our approximate treatment of the interaction (and 
also in bosonization) are taken into account. They set a new upper energy scale B* 
for the simple scaling behavior considered here which is presumably smaller than the 
band width B. Power laws might further be obscured by the presence of relaxation and 
dissipation in the leads, not included in our calculations. 

The optimal setup to observe the boundary exponent aend is an end-contacted, 
long wire with no overhanging parts of wire or leads for which one can expect scaling 
to hold for vp/N <^ T <^ B*. Although our approximation does not include the bulk 
LL exponent the results indicate that the best setup to observe ctbuik is a long wire with 
leads which both couple to the bulk of the wire and which terminate at the contacts. 

Our approach allows for more general systems with e.g. leads of asymmetric widths, 
random wire-lead couplings or screened interactions in the contact regions. 
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Appendix 

Calculation of resolvent matrix elements 

In order to compute the right-hand side of the flow equations f|T8|) and f|T9|) for the 
self-energy, one needs the tridiagonal elements of the Green function Q. In the case 
of Id leads itself is tridiagonal and these elements can be computed in order N 
time fl6[ [28] with being the length of the wire. In the geometry with stripe-like 
leads, introduced in section 13. ![ the inverse Green function is of the structure shown in 
figure lAli Large tridiagonal blocks (the middle and overhanging parts of the wire) are 
connected on the sub- and superdiagonal to (small) full blocks (the contact regions). 
We can exploit the tridiagonal structure of large parts of this matrix enabling us to 
treat lattice systems with up to 10^ site^. To achieve this, we make multiple use of the 
projection formalism introduced in section 13.11 to split the problem such that we only 
have to invert matrices of size and structure of the isolated blocks 1 to 5 (figure RTl) . 
This allows us to use the effective order A^ time algorithm for the tridiagonal parts and 
some standard algorithm [27] for the (small) full matrices. Thus, the bottleneck of this 
algorithm is the size of the full matrices. 

I Our calculations were performed on a standard desktop PC. Thus the maximum system size reachable 
can be easily extended by using more elaborate hardware. 
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Figure Al. Structure of the inverse Green function for stripe-like leads: block 1,3 and 
5 are tridiagonal (left, right overhang and middle part of wire), block 2 and 4 are full 
matrices (left and right contact region), the dots represent the connecting elements; 

{Nj , Nj ) is the lower right element of the j-th block 



In the first step of tlie algoritfim tlie full Hamiltonian h appearing in the Green 
function Q{z) = [z — h)^^ is split into five smaller parts hj, referring to block j G 
{1, ..,5} respectively, and elements hi'^^^ = (3jj^i\Nj){Nj + 1| + H.c. with j G {1, ..,4} 
connecting block j and j + 1 (with G M), i.e. 

5 4 

h = Y.^, + Y.^^ot (A.l) 

i=i i=i 

We next define the projection operators Pj on each block. Then we project the full 
Green function Q onto the block, whose elements we want to compute, thereby reducing 
the effort of inverting the full matrix to the inversion of an effective matrix of size and 
structure of block j of the inverse Green function. We will exemplify this by calculating 
the elements of block 1 as well as those elements on the sub- and superdiagonal 
connecting block 1 and 2. Making subsequent use of the projection formula (ITSl) . we 
obtain 



PiQPi = 

(z - hi 

P2Q^'^P2 = 

[z - h2 

(^z - h-s 

PaQ^^^Pa = 
{^z — hi 

P.Q^^'^P, = 



(^h\N2){N2 + l\g^^^\N2 + 1){N2\ 
z - h^Y^ 



-1 
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with Q^^^ = (^z — Yl^j=i + Si=i Hon^^ ■ We calculate the elements of block 1 of the 
Green function by going through this hierarchy from bottom to top, having only to invert 
matrices of size and structure of the isolated blocks at each step. The computations for 
the other blocks of Q can be performed similarly. 

The element {Ni + ll^lA?"!) can be calculated by splitting h into h = Yl^j=i + 
Ei=2 hi^n^ and hl^^ to write 

g = g + gh'jj + GhiilGhli^g 

with Q = (^z — h^ . After projecting Q onto the relevant block P1QP2, one ends up 
with 

. + 11^1^ ^ = MN^\g\N^){N, + l\g\N, + l) 

X (l+/322(iVi|^|iVi)(iVi + l|^|iVi + l)) . 

The relevant elements of g can be efficiently calculated using the algorithm introduced 
above for the blocks. Again the calculation of the other elements {Nj + ll^lA'j) can be 
performed analogously. 

For the calculation of the transmission probability (fT2l) we need all matrix elements 
{f\g\j) with f e {Ni + 1, .., N2} and j G {N5 + 1, .., N4} of the Green function g-\ We 
use an algorithm which only requires the inversion of matrices of size and structure of 
the isolated blocks 1 to 5. We split the full Hamiltonian ( lA.ll) into blocks h = Yl^j=i 
and connecting elements /icon = X]j=i Hon^ decompose 

g = g + ghcong + ghconghcong 

with g = (^z — . Multiplication with projection operators P2 and P4 yields 

P^gP, = P,g (/?2,3|iV2)(A^2 + l\+f3l,2\Nl + I) {Nil) 

g (/53,4|iV3)(iV3 + 1| + A,5|iV4 + l)(iV4|) ^^2 • 

The elements of g can be easily computed as the disconnected blocks of g^^ can 
be inverted individually. The elements {N2 + l\g\N3), (A^il^jiVa), {N3 + l\g\N4), 
(A''i|^|A'4 + 1) can be calculated by splitting h into blocks and bonds and projecting onto 
the blocks needed. Details of the rather lengthy procedure can be found elsewhere [22] . 
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